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> ' Abstract 

in 

We provide analytic proofs for the shape invariance of the recently discovered 
(Odake and Sasaki, Phys. Lett. B679 (2009) 414-417) two families of infinitely many 
exactly solvable one-dimensional quantum mechanical potentials. These potentials 
are obtained by deforming the well-known radial oscillator potential or the Darboux- 
Q^ | Pdschl- Teller potential by a degree I {I = 1,2,...) eigenpolynomial. The shape in- 

variance conditions are attributed to new polynomial identities of degree 3£ involving 
cubic products of the Laguerre or Jacobi polynomials. These identities are proved 
^ , elementarily by combining simple identities. 

a: 

1 Introduction 



in 



In a previous Letter [Tjj, two sets of infinitely many shape invariant [2] potentials were de- 
rived by deforming the radial oscillator potential [31 H] and the Darboux-Poschl- Teller (DPT) 
potential [5j [6] in terms of a polynomial eigenfunction of degree £(£ = 1,2,...). As the main 
part of the eigenfunctions of these exactly solvable quantum mechanical systems, the excep- 
tional (Xe) Laguerre and Jacobi polynomials were obtained [T|. The lowest (I = 1) examples, 
the X\ Laguerre and Jacobi polynomials, are equivalent to those introduced in the pioneer- 
ing work of Gomez-Ullate et al [3, [8] within the Sturm-Liouville theory. The reformulation 
in the framework of quantum mechanics and shape-invariant potentials was done by Quesne 



et al [9j [10]. By construction these new orthogonal polynomials satisfy a second order dif- 
ferential equation (the Schrodinger equation) without contradicting Bochner's theorem 
since they start at degree £, (£ = 1,2,.. .) instead of the degree zero constant term. 

Here we present analytical proofs of the main assertion that these deformed potentials are 
indeed shape invariant, which could not be given in the Letter [1] due to the lack of space. 
For both the deformed radial oscillator potential and deformed trigonometric/hyperbolic 
DPT potentials, the condition for shape invariance for each I is satisfied if a certain identity 
involving cubic products of the Laguerre or Jacobi polynomials holds. To the best of our 
knowledge, these infinitely many identities were not presented before. We show that these 
identities are derived by combining several elementary relations among the Laguerre or Jacobi 
polynomials. 

This paper is organised as follows. In section two the general setting of shape invariance 
and the deformation of a potential in terms of a polynomial eigenfunction is recapitulated 
together with many preparatory materials. In section three we show that the condition 
for shape invariance for each I is attributed to a certain polynomial identity of degree 3£ 
involving cubic products of the Laguerre or Jacobi polynomials of various parameters, (13.41) 
and (13. lip . Then we show that these identities are simply derived by combining simple 
relations among the Laguerre polynomials (I3.5p - (I3.6I) and Jacobi polynomials (I3.12p - (l3.13p . 
The final section is for a summary and comments on related results to be published in a 
forthcoming paper. 

2 General setting 

As in [JJ, let us start with a generic one-dimensional quantum mechanical system in a 
factorised form: 

H = AA = p l + U(x), p = -id x , U(x) = (d x w(x)) 2 + d 2 x w(x), (2.1) 
A = d x -d x w(x), A = -d x -d x w(x). (2.2) 

Here we call a real and smooth function w(x) a prepotential . It parametrises the groundstate 
wavefunction 4>q(x), which has no node and can be chosen real and positive, (j)o(x) = e w ^ x \ 
It is trivial to verify A<po(x) = and Tiip^x) = 0. 

Shape invariance is realised by specific dependence of the potential, or the prepotential 
on a set of parameters A = (Ai, A2, • • •), to be denoted by w(x; A), -4(A), H(X), S n (X), etc. 
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The shape invariance condition to be discussed in this paper is written simply: 

^(A)^(A) t = A(X + <5)U(A + 6) + £ X (A), (2.3) 
or (9^(x;A)) 2 -^«;(x;A) = (d x w(x;\ + 6)) 2 + d£w(x;\ + 6) +S 1 (\), (2.4) 

in which 8 is a certain shift of the parameters. Then the entire set of discrete eigenvalues 
and the corresponding eigenf unctions of H = TC(X), 

H(\)<p n (x- A) = £ n (\)Mx; A), n = 0, 1, 2, . . . , (2.5) 

0n(^; A) = (f) (x; X)P n (r)(x); A), (2.6) 

is determined algebraically [21 [U [121 [13] : 

n-l 

£ n (A) = ^2e 1 (\ + kS), (2.7) 

fc=0 

n (x; A) oc ^(A)U(A + <5)t ■ ■ ■ A(X + (n - l)<5) f x e ™(*;A+n<s). (2.8) 

The polynomial eigenfunction P n (r](x);A), which is the Laguerre or Jacobi polynomial in 
T](x), satisfies 

- d 2 x P n ( V (x); A) - 2d x w{x; X)d x P n {rj{x); A) = £ n (X)P n (rj(x); A). (2.9) 



Here rj(x) is a function of x called the sinusoidal coordinate 

In pQ a shape invariant prepotential w(x;X) = Wo(x; A) is deformed by a polynomial 
eigenfunction ^ of its Hamiltonian to produce another shape invariant prepotential we (£ = 
1,2 



f \^ def \j ex\ i i 6(^(a;); A + <5) 

A) = w (a:; A + id) + log ...... , (2.10) 

in which ^ is related to the polynomial eigenfunction P n above (12. 6p . It should be noted 
that the normalisation of the polynomial A) is irrelevant to the deformation. The I = 
case corresponds to the original system. The £-th Hamiltonian and eigenfunctions, etc are 
given by 

At(X) = d x - d x Wi{x; A), A(A) f = -d x - d x w e (x; A), (2.11) 
H e (X) d = Ae(X)U t (X), (2.12) 
H e {X)(f> t>n {x; A) = £ n {X + I6)<j>t, n (x\ A), (2.13) 



3 



def e w o( x '> x + tS ) 

4>e, n {x; A) = tp e (x; A) P e , n (rj(x); A), if> t (x; A) = — (2.14) 

The orthogonality of the eigenfunctions of the Hamiltonian 7ie(X) reads 

y ^(s; A) 2 P^ n (r/(x); X)P £>m (T)(x); X)dx = h^ n (X)5 nm , h t , n (X) > 0, (2.15) 

in which ^(x; A) 2 is the orthogonality measure and Pg >n (x; A) is the n-th member of the 
exceptional (Xg) orthogonal polynomial. It is expressed in terms of ^(x)'s and P n (x)'s as 
shown in f[2T2"2"j) . (12T2TD and fl2TMj) . 

In this paper we will demonstrate that the deformed prepotential wg (£=1,2,.. .) actually 
satisfies the shape invariance condition 

A,(x;A) = 0, (2.16) 

in which Ai(x; A) is defined by 

A e (x;\) = (d xWi (x; X)) 2 -d 2 x w e (x; \)-(d x w t {x; X+8)) 2 -d 2 x wg{x- X+8)-S 1 (\+£8), (2.17) 

for the three cases, the radial oscillator and the trigonometric/hyperbolic DPT discussed in 
PP. The proof consists of two steps. Firstly in section three we transform, by utilising the 
differential equations for £i(r)(x)', A) (12.411) etc, the shape invariance condition (12.171) into an 
identity involving products of three Laguerre or Jacobi polynomials of various parameters, 
(13. 4p and (13.111) . The trigonometric and hyperbolic DPT lead to the same identities. Sec- 
ondly these cubic identities are proved by combining simple identities among the Laguerre 
or Jacobi polynomials of neighbouring degrees n, n — 1 and parameters a, a ± 1, (3, (3 ± 1. 
Here we show various data necessary for the proof. They are recapitulated from pp. 



1 radial oscillator: 

X = g, 6 = 1, g>0 1 (2.18) 

£ n (X) = An, n{x) = x 2 , < x < oo, (2.19) 

(t) {x; A) = e~~x 9 w (x; A) = -y + g log a;, (2.20) 

P n (x;X)= i L { r l2 \x), 6(*;A) = f 4 9+ ^ } (-*), (2-21) 

P iin (x; A) = g + l)P n (x; g + i) - ^{x; g + 2)P n _ 1 (x; g + £). (2.22) 
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2 trigonometric DPT: 



X=(g,h), S =(1,1), h>g>0, (2.23) 

£„(A) = 4n(n + g + ft), rj(x) = cos 2x, < x < ~, (2.24) 

0o(x; A) = f (sin a;) 9 (cos x)' 1 <^> iyo(a;; A) = g log sin a; + ft log cos x, (2.25) 

P n (x; A) ^P B °'-*' h -* ) (z), &(s;A) d = f Pi'^'^V), (2.26) 

Pt, n (x) A) = a*, n (x; A)P n (x; X + 16) + b t , n {x] X)P n ^(x; X + £S), (2.27) 

/ \\dfif t / Ll , , ,x , 2n(-# + ft + £ - 1) &_i(a:; g, ft + 2) 
a^ n (x; A) = 5- + 1, ft + 1) + 



(-g + ft + 2£ - 2)(g + ft + 2n + 2£ - 1) 
n(2ft + 41-3) &_ 2 (x; g + 1, ft + 3) 
(2# + 2n+l)(-# + ft + 2£-2) ' 



3 hyperbolic DPT: 



(-g - ft + 2£ - 2)(g - ft + 2n + 2£ - 1) 
n(-2ft + 4£ - 3) &- 2 (z; g + 1, ft - 3) 
(2# + 2n+ !)(-#- ft + 2£- 2) ' 



(2.28) 



7 , ^def (-<7 + h + l-l)(2<7 + 2n + 2l-l) 

MX; A) = (2, + 2n + l)(, + ft + 2n + 2£-l) * * + 2) ' (2 ' 29) 



A= f (s,ft), S = (1,-1), ft>^>0, £<n B = f [§(ft -<?)]', (2.30) 

£ n (A) = 4n(ft — g — n), rj(x) = cosh2x, < x < 00, (2.31) 

4>o(x; X) ^ (sinhx) g (cosh.x)~ h ^ wq(x; X) = g log sinh x — ft log cosh x, (2.32) 

P n (x;X) d ^Pt h ' h ' l2 \x), &(x;X) = f P^-*-™'^), (2.33) 

P e>n (x; A) = a/, n (x; A)P n (x; X + £S) + b e>n {x; A)P n _!(x; A + IS), (2.34) 

/- ^ def t^ 1 1 a ^1 2n(-^-ft + £-l)6_i(x;^,ft-2) 
a £jn (x; A) = g- + 1, ft - 1) + 



(2.35) 



7 , . x def ( -g - ft + £ - 1 ) (2g + 2n + 21 - 1 ) 

A) = (2g + 2n + l)(g-fr + 2n + 2l--l) * " ^ ^ 

In ( 12 . 30 j) [x\ denotes the greatest integer not equal or exceeding x. Here Ln*\x) is the 
Laguerre polynomial and Pn (x) is the Jacobi polynomial. The polynomial eigenfunction 
Pe, n {x; A) (I2.22p . (12.271) or (12.341) is a degree £ + n polynomial in x. In a future publication 
[T5] we will present equivalent forms of the polynomials P^ n (a;; A) which appear much 
simpler than those given in (12.221) . (12.271) and (12.341) . Needless to say that the deforming 
polynomials £e(T](x); A) and £4(77 (a:); A + <5) are of the same sign in the domain for all the 
three cases. In other words, the deforming polynomials do not have a zero in the domains 
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listed in (I2.19p . (I2.24p and (12.311) . respectively. To see this we use the expansion formula of 
the Laguerre and Jacobi polynomials: 



x 



-n k 



a + k + l)n- k x h , 



1 n 

-T 

n\ ^ k\ 

k=0 

(a + l)„>A 1 (-n) k (n + a + (3 + l) k fl-x\k 



k=0 



(2.37) 
(2.38) 



For the radial oscillator we obtain 



e*,(*);A) = £ to+ .'+ t T.? ) '- t *»>0. 



k=0 



k\ (£ — k)\ 



(2.39) 



and for the trigonometric/hyperbolic DPT 



{ ^ + ^^^-^^-^-^(sinxr>0, 



k=0 



k\ (g + e-k+i) t 



(2.40) 



k=0 



k\ (g + £-k + \) k 



because each term in summation is positive. This guarantees the positive defmiteness of 
the orthogonality measure ipg(x; A) 2 (12.141) and the singularity free structure as well as 
the hermiticity (self-adjointness) of the Hamiltonian. The oscillation theorem for the one- 
dimensional quantum mechanical systems dictates that the n-th excited state polynomial 
eigenfunction Pi jn (r](x); A) has n zeros in the domain (12.191) . (12.241) or (12.311) . although it is 
a degree £ + n polynomial in rj. 

Let us note that the terms of the form of the complete square (d x £e(ri(x); A) / £i(r](x); A)) 2 , 
{d x ^(r](x);X + 6)/^(r](x);X + 6)) 2 and (d x £ £ (r)(x)\ A+ 2d)/&(r)(x); A + 2<5)) 2 , cancel out in 
Ae(x; A) (12. 1 7D . Then we use the fact that, corresponding to (\2.9L the deforming polynomial 
£e(rj(x); A) also satisfies a second order linear differential equation 



d l&(v(x); A) ~ ^d x w (x; A, £)d x £ e (ri(x);\) = £ e (X)^(r](x); A), 



(2.41) 



where w (x; \,£) and £e(\) are 



w (x] A,. 



1™2 
2 



(g + £ — 1) logx 
w (x; —g — £, h + £ — 1) 
wq(x; —g — £,h — £ + l) 



radial osci. 
trig.DPT 
hyper. DPT, 



(2.42) 
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{— 4£ : radial osci. 

£ e (-g-£,h + £-l) : trig.DPT (2.43) 
£ t (-g-e,h-e+l) :hyper.DPT. 

Substituting (12.101) into (I2.17P and using the shape invariance of the undeformed system 
A (a;; A + 18) = and (j2"H]) . we obtain 

A £ (x; X) x ~ Ci(v, A) A + S) UV, A + 28) 
= d x w (x; X + 28, I) d x Uv, A + 26) Uv, A) Uv, X + 8) 

- d x w (x; X, £) dg-Ceim A) Zt{m A + <5) A + 28) 

+ \ (£ e (X + 28) - £ t (X)) Uv, A) Uv, A + 8) Uv, A + 28) 

+ d x w (x; X + £8) (d x Uv, A + 8) Uv, A) - d x Uv, A) Uv, A + 8))Uv, A + 25) 

- d x w (x; X + £8 + 8) (dMv, A + 28)Si{w X + 8)- d x &{v, A + %(tj; X + 28))&{ m X) 

- dMv, A) d x Um a + s) Uv; a + 28) + d x Um a + s) dMv; a + 28) Urn A), (2.44) 

where rj = i](x). It is essential that the r.h.s. of (12.441) is now a polynomial in f](x), because 
d x rj(x)d x WQ and d x r](x)d x w and (d x r](x)) 2 are expressed by rj(x), 







' g-v( x ) 


radial osci. 




d x r](x)d x w {x; X) 


= 2 x < 


~{g-h + (g + h)r](x)) 


trig.DPT 


(2.45) 






K g + h + (g - h)r](x) 


hyper. DPT, 








' g + i-\+r\{x) 


radial osci. 




d x T](x)d x w (x; X,£) 


= 2 x < 


g + h + 2£-l + (g-h + l)r)(x) 


trig.DPT 


(2.46) 






k -(g-h + 2£-l + (g + h + l) V {x)) 


hyper. DPT, 




(d x 7](x)) 2 




' T](x) 


radial osci. 




= 4 x < 


l-^xf 


trig.DPT 


(2.47) 








hyper.DPT. 





3 Cubic identities 

The conditions for the shape invariance (I2.16P are shown to be equivalent to cubic identities 
involving the Laguerre or Jacobi polynomials. 

3.1 radial oscillator 

We fix £ and use a new parameter a instead of g, 



By using the forward shift relation for the Laguerre polynomial, 

d x L^(x) = -Lt + 1 1 \x), (3.2) 

the polynomial ^ and its derivative are expressed as 

UV, A) = ^(-rj), d x Uv, A) = d xV iS&i-rj), (3.3) 

where 77 = r](x). After replacing —r](x) with x and dividing by 4, the condition for the shape 
invariance (I2.16P with (12.441) . is transformed into the following polynomial identity of degree 
3£ in x which contains products of three Laguerre polynomials of various parameters: 

= - xL { tX 2 \x)L { t~ l \x)Lf\x) - aL { t\(x)Lf +1 \x)Lf\x) 
+ (x + a + l)4 a + 1 \x)L[ a+1 \x)4 a - 1 \x) 

+ xLf\{x)Lf^\x)Lf +1 \x) - xLtY\x)L ! f_ + 1 2) (x)L ! f- 1 \x). (3.4) 

To the best of our knowledge, this identity has not been reported before. For £ = this 
identity is trivial, since L^\{x) = 0. For lower I it can be easily verified by direct calculation. 

Below we will prove the identity ( 13. 41) for an arbitrary positive integer i by combining 
a few elementary relations among the Laguerre polynomials of neighbouring degrees n and 
n — 1 and neighbouring parameters a, a ± 1: 

Lemma (A) ^(x) + L^x) = L^(x), (3.5) 

Lemma (B) xL { ^\x) - aL { °\(x) = -n^ a_1) (x), (3.6) 

which can be verified elementarily based on the expansion formula of the Laguerre polynomial 
(I2.37p . The r.h.s. of the identity (13. 4p can be written as 

r.h.s. of (J32D = -xLt\ 2 \x)Lf- 1 \x){Lf\x) + Lt + l 1 \x)} 
+ xLl^^Lf^ix^Lt^ix) + Lt\{x)} 
+ L[ a+1 \x)((a + l)^ 1 ^)^" 1 ^) - aLt\{x)Lf{x)) 
= Lf +X \x)( K Lf\x){xLf^\x) - aLf\{x)} 

+ Lf~ l \x){-xLtf\x) + (a + 1)^%)}) 
= L^ 1 \x)(L^\x)(-iLt 1 \x)) + Lf- l \x)(+lLt\x))) = 0. (3.7) 

Here we have used Lemma (A) in the first two curly brackets {■ ■ ■ } and Lemma (B) in the 
next two curly brackets. This concludes the proof of the identity (13. 4p . 
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3.2 trigonometric DPT 



We fix t and use new parameters a and (3 instead of g and h, 



def 

a = 



-g-£-\, p = h + t-\. 



(3.8) 



By using the forward shift relation for the Jacobi polynomial, 




(3.9) 



the polynomial ^ and its derivative are expressed as 



A) = Pp P) (v), dJU{m \) = d m ri\{t + a + P + l)Ptr P+L \vl (3-10) 



where 77 = r](x). The condition for the shape invariance (12.161) with (12.441) . after replacing 
rj(x) with x and dividing by —(£ + a + (3 + 1), is simplified to a polynomial identity of degree 
3£, which contains products of three Jacobi polynomials of various parameters: 



To the best of our knowledge, this identity has not been reported before, either. For £ = 
this identity is trivial, since P*f{^\x) = 0. It is straightforward to verify this identity for 
lower I by direct calculation. 

Below we will prove the identity (13.111) for an arbitrary positive integer I by combining a 
few elementary relations among the Jacobi polynomials of neighbouring degrees n and n — 1 
and neighbouring parameters a ± 1, /3 ± 1: 



2(a - l)P { n a ~ m {x) -{n + a + (3){l - x)P^ +1 \x) = 2(n + a - l)P^- 2 ^ +1 )(x), (3.12) 



2{f3 + l)Pt- 1 ' ni \x) + (n + a + (3 + + x)Pi a J +2 \x) =2(n + (3 + l)P^\x). (3.13) 




(3.11) 



Lemma (C) 



Lemma (D) 
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It is straightforward to demonstrate Lemma (C) by using the expansion formula for the 
Jacobi polynomials (j2.38|) . By using the property Pn®(— x) = (—l) n Pn ,a \x), Lemma (D) 
is obtained from Lemma (C) with the replacements x — > —x, a — >/3 + 2, (5 — > a — 1. The 
r.h.s. of the identity (13. lip can be written as 

r.h.s. of (EUD = (1 + x)p£~ 1 ' 0+3 \x)P^\x) 

x {2(a - l)P^ +1 \x) - {£ + a + /3 + 1)(1 - x)P^f +2 \x)} 

+ (1 - x)P£l^ +1 \x)Pt 2 '" +2 \x) 

x{2((3 + l)P^ p+l \x) + {£ + a + (3 + 1)(1 + x)P^f +2 \x)} 



2(a(l + x) + (J3 + 2)(1 - x))p£f +2 \x)P^\x)Pl a - 2 ^ +2) ~ 



x 



= P^\x)Pt 2 '" +2 \x) 

x((l+x){2(£ + a- l)p£- 1JJ+3) (a;) - 2ai£f +2) (x)} 

+ (1 - x){2(£ + /3 + l)P^+ 1,/3+1) (x) - 2(/3 + 2)P, ( !f +2) (x)}' 

x (£ + a + p + 2)((1 + - 1) + (1 - x){x + 1)) 
= 0. (3.14) 

Here we have used Lemma (C),(D) in the curly brackets {•■•}■ This concludes the proof of 
the identity (pTTIjl . 

It is well known that the Laguerre polynomial is obtained from the Jacobi polynomial in 
the following limit: 

hmPt^{l-f)=Li a \x). (3.15) 

When the same limit is applied, the Lemma (C)-(D) reduce to Lemma (A)-(B). Likewise 
the cubic identity for the shape invariance of Xi Jacobi polynomial (13.111) reduces to that of 
the Xi Laguerre polynomial (13 ,4p . when divided by —4 and a is replaced by a + 1. 

3.3 hyperbolic DPT 

Here we briefly remark that the shape invariance condition for the deformed hyperbolic DPT 
reduces, as expected, to the same identity as that for the trigonometric DPT ( 13. lip derived 
above. We fix I and use new parameters a and f3 instead of g and h, 

a d ^-g-£-l, p^-h + i-l (3.16) 
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By using (13.91) . the polynomial & and its derivative are expressed just same as (13.101) . The 
conditions for the shape invariance (12.161) with (12.441) . after replacing r](x) with x and dividing 
by £ + a + (3+ 1, is simplified to the same identity (13. lip as that for the trigonometric DPT. 

4 Summary and Comments 

Analytic proofs are provided for the shape invariance of the recently derived infinite family 
of potentials [lj obtained by deforming the radial oscillator potential and the trigonomet- 
ric/hyperbolic Darboux-Poschl- Teller potential [3J E] by a degree £ (£ = 1,2, . . .) eigenpoly- 
nomial. The shape invariance conditions are reduced to new polynomial identities of degree 
3£ involving cubic products of the Laguerre (13 .4p or Jacobi polynomials (13. lip . Then these 
identities are proved elementarily by combining simple linear identities (I3.5l) - (l3.6p among the 
Laguerre and (I3.12I) - (I3.13I) among the Jacobi, polynomials of neighbouring degrees n, n — 1 
and of neighbouring parameters, a, a ± 1 and /3, /3± 1. Even these linear identities seem not 
widely recognised. 

The totality of the eigenvalues and the corresponding eigenfunctions together with the 
normalisation constants etc. of these infinite family of quantum mechanical systems are 
obtained via the Rodrigues type formulas (12.71) and (12. 8p and reported as (17)-(21), (29)-(36) 
and (44)-(46) of [I]. In a future publication [13] ; we will present analytic derivation of various 
results reported in pQ. They include: derivation of equivalent but much simpler looking forms 
of the Xi polynomials together with the normalisation constants, the verification of the 
actions of the forward and backward shift operators on the X% polynomials, Gram-Schmidt 
orthonormalisation for the algebraic construction of the Xi polynomials, the analysis of 
the second order differential equations for the Xi polynomials within the framework of the 
Fuchsian differential equations in the entire complex x-plane, etc. In these analysis, the linear 
identities Lemma (A)-(D) play important roles. The forward and backward shift relations 
mentioned above (eqs.(49),(22),(37),(47) of [1]) are also reduced to polynomial identities 
involving cubic products of the Laguerre or Jacobi polynomials. We will provide proofs for 
them in the future publication [13] . too. 

Let us also mention that the same method, deformation in terms of a degree £ eigen- 
polynomial, applied to the discrete quantum mechanical Hamiltonians for the Wilson and 
Askey- Wilson polynomials produced two sets of infinitely many shape invariant systems 
together with exceptional (Xi) Wilson and Askey- Wilson polynomials (£ = 1,2, . . .) [T6] . 
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